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Simultaneous  photographs  are  taken,  of  the  same  aerial  point,  by- 
two  cameras  on  the  ground.  Stars  also  are  photographed,  on  some  or  all 
of  the  plates.  Here  the  mathematical  procedure  is  developed  and  described 
for  finding  accurately  the  position  of  the  aerial  point  with  respect  to 
a  set  of  terrestrial  axes,  g  ,  n  ,  end  £  .  The  formulae  are  designed 
to  expedite  the  numerical  calculations  as  much  as  possible,  consistently 
with  an  accuracy  of  the  order  of  a  second  of  arcj  i.e.,  of  better  than  a 
foot  at  50,000  feet.  Another  less  accurate  procedure  is  also  described 
that  does  not  involve  star-images  on  all  of  the  plates.  The  present  dis¬ 
cussion  is  applicable  to  the  case  where  both  camera  axes  are  vertical. 

Some  of  the  formulae  are  applicable  also  to  oblique  camera  axes. 


1.  Introduction.  In  order  to  obtain  the  fundamental 
data  from  which  bomb  ballistic  tables  can  be  prepared,  it 
is  necessary  to  know  the  position  and  velocity  of  an  airplane 
at  the  instant  when  it  releases  a  bomb.  Both  the  position 
and  velocity  can  be  determined  from  a  knowledge  of  the  posi¬ 
tion  of  the  airplane  at  each  of  several  known  times  and  from 
a  knowledge  of  the  time  of  release.  Nearly  instantaneous 
light  signals  can  be  emitted  from  the  airplane  at  known 
times,  and  each  of  these  signals  can  be  photographed  by  two 
cameras  mounted  rigidly  on  terrestrial  piers.  The  problem 
with  which  we  are  here  concerned  is  the  determination  of  the 
position  in  space  of  each  of  these  light  signals,  from  the 


The  position  of  an  aerial  point  is  determined  when 
its  directions  from  both  cameras  are  known.  In  order  to  be 
of  use  in  range  bombings  from  great  altitudes,  these  direc¬ 
tions  must  be  found  to  a  very  high  accuracy.  The  most  accu¬ 
rate,  and  yet  not  prohibitively  laborious  or  lengthy  pro¬ 
cedure  is  to  photograph  known  ’’comparison"  stars  on  the  same 
plates  that  photograph  the  aerial  points.  Then,  by  adapting 
to  the  present  problem  the  simple,  mechanical,  and  yet  ex¬ 
ceedingly  accurate  methods  of  photographic  astrometry,  one 
can  find  the  directions  of  the  rays  to  the  aerial  points 
from  a  knowledge  of  the-  positions  of  the  comparison  stars, 
accurately  listed  in  star  catalogues.  The  x  and  y  coordi¬ 
nates,  of  the  images  of  the  stars  and  of  the  aerial  points 
on  the  photographic  plates,  must  be  measured  on  a  measuring 
engine.  The  star-images  determine  a  set'  of  ’’plate  c oils t ants", 
and  the  plate  constants  in  turn  determine  the  direction 
ratios  to  the  aerial  points.  It  is  not  necessary  for  the 
ways  of  the  engine  to  be  exactly  perpendicular  to  each  other 
or  parallel  to  the  edges  of  the  plate;  it  is  not  necessary 
for  the  scales  of  the  measuring  engine,  in  x  and  y,  to  be  the 
same;  no  fiducial  marks  are  necessary  on  the  plates;  the 
focal  lengths  of  the  cameras  need  not  be  accurately  known; 
and  errors  in  the  assumed  orientations  of  the  cameras  pro¬ 
duce  only  second-order  errors  in  the  resulting  directions 
to  the  aerial  points. 

An  alternative  and  less  accurate  procedure  is  to  / 
determine  the  plate  constants  on  some  plates  from  the  star- 
images,  and  to  find  the  directions  of  the  rays  to  an  aerial 
point  on  other  plates  by  employing  fiducial  marks  impressed 
upon  the  plates  by  the  camera.  The  alternative  procedure 
obviates  the  need  for  employing  star-images  on  all  the 
plates,  but  presumes  a  very  high  degree  of  fixity  and  per¬ 
manence  in  the  mountings  and  orientations  of  .the  cameras. 

Here  we  develop  and  present  the  necessary  formulae, 
many  of  which  are  mere  adaptations  of  the  well-known  formulae 
of  photographic  astrometry.  Astrometry,  however,  concerns 
the  relative  directions  of  stars  with  respect  to  each  other; 
in  the  present  problem  we  have  also  to  find  directions  rela¬ 
tive  to  the  surface  of  the  earth.  The  present  problem,  there¬ 
fore,  has  some  new  aspects.  Further,  it  is  desirable  to 
simplify  some  of  the  conventional  astrometric  formulae  with  a 
view  to  accelerating  routine  computation,  to  diminishing  the 
frequency  of  errors  of  computation,  and  to  facilitating  the 
checking  of  the  various  stages  of  the  reductions.  To  avoid 
errors  arising  from  the  entering  of  tables  incorrectly,  an 
effort  will  be' made  to  avoid  the  use  of  trigonometric  func¬ 
tions,  by  using  direction  cosines  instead.  After  some  pre¬ 
liminary  transformations  have  been  made  to  direction  cosines, 
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no  further  use  will  be  made  of  trigonometric  functions.  It 
is  believed  that  the  new  forms  into  which  some  of  the  old 
astrometric  formulae  have  been  recast  may  be  useful  outside 
of  the  present  problem  in  astrometry. 

The  use  of  two  terrestrial  cameras  was  urged  by 
Professor  Henry  Norris  Russell,  to  whom  the  writer  is  grate¬ 
ful  also  for  his  valuable  and  specific  suggestions.  Among 
them  were  the  finding  of  lens-distortion  by  examining  the 
residuals,  of  a  least  squares  solution  for  plate  constants. 

He  also  pointed  out  that  whenever  the  "base"  of  a  plate 
is  the  zenith,  then  the  direction  cosines  in  the  altitude- 
azimuth  system,  to  any  point,  can  be  found  at  once  from 
its  standard  coordinates. 

We  take  up  the  problem  in  the  order  in  which  the 
computations  should  be  made.  We  discuss  first  the  reduction 
of  mean  star  positions,  found  from  a  star-catalogue,  to 
apparent  positions.  Here  the  most  expeditious  procedure, 
since  standard  coordinates  must  be  computed,  is  to  compute 
direction  cosines  at  once  and  then  correct  the  direction 
cosines,  rather  than  the  right  ascensions  and  declinations, 
for  precession  since  the  beginning  of  the  year,  for  nutation, 
and  for  aberration.  Then  standard  coordinates  are  to  be 
computed,  to  be  corrected  for  lens-distortion  and  third-order, 
terms  in  refraction.  Then  the  plate  constants  are  found,  ana 
from  them  and  the  plate-coordinates  of  the  aerial  point  one 
computes  the  direction  ratios  of  the  rays  to  the  aerial 
point.  From  these  (corrected  for  refraction  and  distortion),., 
the  terrestrial  coordinates  of  the  aerial  point  follow  at 
once. 

In  subsequent  sections  we  show  how  to  follow  the  . 
alternative  procedure  that  does  not  use  star-places  on  all 
the  plates. 

Equations  to  be  used  by  computers  are  numbered  with 
Roman,  other  equations  with  Arabic,  numbers. 


2.  The  Reduction  of  the  Star-Places.  Reduce  the  cata¬ 
logue  mean  positions  of  the  stars  to  the  beginning  of  the 
current  year,  by  applying  the  annual  and  secular  variations 
(which  include  proper  motion)  listed  in  the  catalogue.  Apply 
also  the  proper  motion,  if  it  amounts  to  more  than  0”.l  (or 
0s. 01)  from  the  beginning  of  the  year  to  date,  and  apply 
further  the  small  diurnal  aberration.  The  equations"  are 
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a  =  a  + 
'  0 


4  0s, P2 


where 


At.2£  tut' 

200 

6=6  +  6t  4  4  u't' 

0  200 


(I) 


a,  6 


a  .6 
o’  0 

*  A 

a,o 

II  11 

a,  6 


t» 


are  the  annual  proper  motions  in  right  ascension 
and  declination, 

are  the  right  ascension  and  declination,  for  the 
mean  equinox  of  the  beginning  of  the  current  year, 
and  epoch  of  date,  affected  by  zenithal  diurnal 
aberration, 

are  the  catalogue  right  ascension  and  declination 
(for  the  mean  equinox  of  the  catalogue) 

are  the  catalogue  annual  variations, 

are  the  secular  variations  in  right  ascension  and 
declination, 

is  the  number  of  whole  years  elapsed  since  the 
epoch  of  the  catalogue, 

is  the  fraction  of  the  current  year  elapsed, 
and, 

0s. 02  is  the  zenithal  diurnal  aberration. 


NOTE:  The  symbols  t  and  ^  will  be  used  later  in  different 
senses. 


It  will  be  noted  that  t  and  t2/200,  as  well  as  t*  ,  are  the 
same  for  all  the  stars  used  and  may  be  computed  once  and  for 
all  for  any  date.  The  above  reductions  are,  therefore,  very 
rapidly  performed  If  t,  t2/200,  and  V  be  entered  on  a  slip 
of  paper.  If  the  Boss  General  Catalogue  is  used,  t.  for  1942 
is,  bf  course,  -8,  and  iyV200  is  +0.32.  The  terms  in  t’  are 
usually  ignorable. 
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A  word  of  explanation  is  required  about  the  diurnal 
aberration.  Its  effect  upon  the  apparent  position  is 


d&  =,  0".S1  cos  <P  cos  T  sec  & 
=  0"  .31  cor  <P  sin  T  sin  6 


where  T  is  the  hour-angle  and  p  the  latitude.  It  can  thus 
never  exceed  0” .31  in  absolute  amount,  in  either  right 
ascension  or  in  declination.  With  vertical  cameras  the 
stars  photographed  will  be  close  to  the  zenith  —  with  the 
Goerz  cameras,  within  20°  of  it.  At  the  zenith 

Hr  ^  ^ 

da  =  0S. 021,  d  6  =  0 . 

It  follows  from  equations  (2)  that  at  Aberdeen  Proving  Ground, 
the  diurnal  aberration  20°  from  the  zenith  can  differ  from 
its  zenithal  value  by  no  more  than  0" .16.  Thus  the  equations 
(I)  allow  correctly  for  diurnal  aberration  to  within  this 
high  accuracy. 

Compute  equatorial  direction  cosines,  L  ,  M  ,  N  , 
for  each  star  by  the  formulae  '  u 

Lq  =  cos6  cosa 

Mq  =  cos6  sina  (ill) 

N  =  sin  6 .  • 
o  i 


It  will  be  noticed  that  these  cosines  relate  to  the  mean 
equator  and  equinox  of  the  beginning  of  the  current  year,  and 
contain  the  diurnal  aberration,  and  the  proper  motion  since 
the  beginning  of  the  year.  These  cosines  will  thus  never 
need  to  be  computed  more  than  once  in  each  year,  except, 

%>  rarely,  in  the  case  of  large  proper  motions. 

v  It  is  next  necessary  to  correct  these  cosines  for 

precession  since  the  beginning  of  the  year,  for  nutation,  and 
for  annual  aberration.  This  correction  is  most  readily 
carried  out  by  the  formulae 
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(IV) 


1 


L  =  Lq  +  AL 
M  =  MQ  +  AM 
N  =  Nq  +  AN 


where  the  corrections  AL,  AM,  and  AN  are  given  by  the  formulae 


P  -  DLq  -  CM0  -  iN0 

10,000  AL  =  PLo  -  FM0  -  ANQ  -  D 
10,000  AM  =  FLq  +  PMq  +  BNQ  +  C 
10,000  AN  =  AL0  -  BMq  +  PNQ  +  i 


(V) 


In  these  equations  the  coefficients  A,  B,  C,  D,  F,  and  _i  de¬ 
pend  only  on  the  date,  and  should  be  tabulated  to  an  accuracy 
of  0.001  for  convenient  reference  in  the  use  of  the  present 
formulae.  They  are  defined  by 

B  =  0.04848  Be 

C  =  0.04848  Ca 

D  =  0.04848 

G  ■ 

F  =  0.72722  (fe  +  f»e) 

i  =  0.04848  ie 

A  =  0.9717  A0  from  1937  through  I960,  and, 

A  =  0.9716  Ae  from  1961  through  1985, 


where  the  symbols  Ae,  Be,  Ce,  De,  fe  ffe  and  ie  denote  the 
quantities  that  are  listed  in  the  American  Eohemeris  and 
-*  Ifeutical  Almanac,  for  each  day  in  the  year,  in  the  section 
4k  dealing  with  star  reductions,  under  the  names  "A",  nB” ,  rC’T, 

^  nD,T,  T'f,T,  r,fMI,  and  "i”  without  the  subscript  "e".  These 

quantities  are  all  in  seconds  of  arc  with  the  exception  of  A  , 
a  pure  number,  and  fe  and  f*e  which  are  in  seconds* of  time.  e 
A,  B,  C,  D,  in  the  notation  of  the  Almanac  are  Besselian  star 
numbers;  f,  f’,  and  i  are  Independent  star  numbers. 
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The  equations  (IV)  and  (V)  have  been  derived  merely 
by  writing  down  the  first  derivatives  of  L,  M,  and  N  with 
respect  to  a  and  6  ,  and  employing  the  formulae  for  da  and 
d6  that  are  listed  on  page  222  of  the  American  Eohemeris  and 
Nautical  Almanac  for  19*38  under  the  title  ’’Independent  Star 
Numbers”.  The  terms  in  the  proper  motion  can  be  dropped  from 
the  latter  formulae  since  proper  motion  has  been  already 
allowed  for.  Formulae  (IV)  and  (V)  follow  after  some  algebraic 
reductions,  which  the  reader  can  readily  verify.  The  second- 
order  terms  in  AL,  AM,  and  AN  are  of  the  order  of  only 
0.0000001,  or  0” .02,  and  have  been  here  Ignored. 

In  using  equations  (V) ,  it  is  sufficient  to  know 
the  coefficients  only  to  the  nearest  0.001  and  to  use  approxi¬ 
mate  values  of  L0,  M0>  and  N0  rounded  off  to  the  nearest 
0.001.  Then  AL,  AM,  and  AN  will  be  accurate  to  the  nearest 
0.0000001,  which  is  sufficient.  It  is  a  good  plan  to  retain 
seven  digits  in  L0>  Mo?  N0?  L,  M,  and  N  and  throughout  the 
subsequent  stages,  except  where  otherwise  specified  in 
following  sections.  We  have  already  pointed  out  that  the 
cosines  L0,  M0,  N0  can  be  computed  once  and  for  all  for  a 
whole  year  of  observation  except  in  the  rare  cases  of  stars 
of  large  proper  motion.  Likewise,  L,  M,  N  will  stay  effec¬ 
tively  constant  on  any  one  night  and  will  vary  only  slightly 
from  one  night  to  the  next. 

The  computation  of  L,  M,  N  can  be  checked  by 
verifying  that  the  sum  of  their  squares  is  unity,  and  this 
check  should  previously  have  been  applied  to  L0,  M0,  and  N0. 


3.  Standard  Coordinates.  Denote  by  t  the  sidereal 
time  of  the  photograph  of  the  star-images.  This  time  will 
differ  for  the  two  cameras  even  if  they  are  exposed  simul¬ 
taneously,  unless  they  lie  on  the  same  meridian  of  longitude; 
denote  by  1,  m,  n  the  direction  cosines  of  a  star  referred 
to  the  axes  of  the  hour-angle  declination  system.  Denote  the 
ordinary  equatorial  system  of  axes  by  x”,  y”,  z”  with  the  x” 
towards  the  equinox,  and  z”  towards  the  north  pole  of  the  sky. 
Denote  the  hour-angle  declination  system  of  axes  by  x’ ,  y’ ,  zf, 
with  x’  pointing  towards  the  intersection  of  the  equator  with 
the  observer’s  meridian,  and  with  the  z’  axis  continuing  to 
point  to  the  north  pole.  The  primed  system  of  axes  is  obtained 
by  rotating  the  unprimed  system  positively  through  the  angler  , 
and  thus 
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(VI) 


\ 

i 


1  =  L  cos  t  +  M  sin  i 
m  =  L  sin  t  +  K  cos  t 
n  =  N 


The  factors  cos  t  and  sin  t  must  be  looked  up,  and  then  all 
the  stars  that  one  is  using  must  have  their  old  cosines, 

L,  M,  and  N  transformed  by  the  equations  (Vi)  into  the  new 
cosines  1,  m,  n.  The  new  cosines  may  be  checked  by  veri¬ 
fying  that  the  sura  of  squares  is  unity,  and  are  thus 
rapidly  computed  and  verified. 

.  It  is  important  (if  one  is  using  the  American 
Eohemeris  and  Nautical  Almanac  for  computing  the  sidereal 
time  from  the  standard  time)  that  the  true  sidereal  time, 
including  terms  of  short  period  in  the  nutation,  must  be 
used  and  not  the  uniform  sidereal  time.  The  "sidereal 
time  of  0^"  listed  in  the  Almanac  may  be  used  properly  for 
the  conversion;  but  not  the  "civil  time  of  sidereal  On",  in 
which  short-period  terms  have  been  ignored. 

By  the  base  of  a  photographic  plate,  we  mean  the 
image  on  the  plate  of  the  infinitely  distant  point  whose 
rays  enter  the  camera  perpendicularly  to  the  plate.  With 
a  horizontal  plate,  the  base  is  the  image  of  the  camera* s 
zenith,  of  which  the  hour-angle  is  zero  and  of.  which  the 
declination  is  equal  to  the  latitude,  <p,  of  the  camera. 

We  are  here  concerned  with  a  horizontal  plate,  whose  base 
is  the  zenith;  in  the  next  section  we  shall  discuss  the 
magnitude  of  the  errors  in  direction  that  can  result  from 
errors  in  levelling,  the  former  errors  being  much  smaller 
than  the  latter. 

Consider  a  set  of  axes  OX,  OY,  OZ  with  OX  parallel 
to  the  y*  axis,  OZ  towards  the  camera’s  zenith,  and  OY  in 
the  plane  ZOz’  and,  therefore,  in  the  terrestrial  north 
direction.  Consider  a  plane  Z  =  K  where  K  is  some  constant. 
The  X  and  Y  coordinates  of  the  intersection  of  this  plane 
(called  the  "standard  coordinate  plane")  with  a  ray  OQ  are 
said  to  be  the  "standard  coordinates"  corresponding  to  the 
direction  OQ.  Let  the  direction  cosines  of  OQ  in  the  x* ,  y* , 
z’  system  be  1,  m,  n.  The  cosines  of  the  angles  between 
the  axes  of  the  two  systems  of  coordinates  are 


X 

Y 

Z 

X* 

0 

-sin  cp 

cos  (p 

y' 

1 

0 

0 

z’ 

0 

cos  (p 

sin  (p 
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I 

Thus  the  direction  cosines  of  OQ  in  the  X,  Y,  Z  system  are 
m,  -  1  sin  <p  +■  n  cos  cp,  1  cos  (p  +  n  sin  <p 
and  the  intersections  desired  are  Z  =  1,  and. 


m 

X  =  K  - r - 

PI  yn 


Y 


=  PJ5 _ i_i 

P  1  +  Y  n 


(VII) 


* where  |3  =  cos  fp,  y  =  sin  cp.  The  factors  p  and  y  are 

mere  constants  dependent  only  on  the  latitude  of  the  camera, 
and  thus  the  equations  (VII)  allow  the  standard  coordinates, 

X  and  Y,  of  any  star  to  be  readily  computed.  For  the  Goerz 
cameras  the  greatest  angular  departure  from  the  OZ  axis  is 
about  20°,  and  it  is,  therefore,  suggested  that  the  constant 
K  be  set  equal  to  10  in  order  not  to  have  too  small  X*  s  and 
Y*s.  Equation  (VII*)  may  be  used  to  check  the  computation 
of  X  and  Y: 

X2  +  Y2  +  K2  =  K2/ (pi  +  yn)2.  (VII*) 

The  equations  (VII)  may  also  be  derived  from  those  of 
Turner*  or  Schlesinger**,  who  employ  a  different  notation,  if 
one  transforms  to  direction  cosines  after  inserting  in  their 
formulae  a  base  whose  declination  equals  the  latitude,  and 
whose  right  ascension  is  equal  to  that  of  the  zenith. 


*  Schlesinger ,  Transactions  of  the  Yale  University  Observatory ,  y  it  Jfo.JO. 
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4*  Accuracy  of  Levelling  the  Plate.  The  standard 
coordinates  will  be  used,  in  a  process  that  is  really  a 
type  of  interpolation,  to  find  the  directions  of  images 
on  the  plate  from  a  knowledge  of  the  directions  of  known 
star-images.  It  is  easy  to  find  approximately  the  effect 
of  small  errors  of  levelling  upon  "the  accuracy  of  direc¬ 
tions  thus  interpolated.  In  the  accompanying  diagram,  C 
is  the  principal  point  of  the  camera  lens,  CD  is  vertical, 
AE  is  the  assumed  horizontal  position  of  the  plate,  BF  is 
its  actual  position,  A  and  E  are  two  star-images  such  that 
AD  =  DE.  Denote  the  error  of  levelling,  namely  the  angle 
BDA,  by  e  ;  denote  the  angle  ACD  by  £  .  Then  A  and  E 


C 


correspond  to  the  standard  coordinates  computed  for  an 
assumed  level  plate.  The  distance  from  B  to  B’  is  nearly 
e  (AD)  where  BT  is  the  foot  of  the  perpendicular  from  B 
to  AE;  and  thus  AB*  is  nearly  e  (AD)  £  .  Now  D  is  the  mid¬ 
point  of  AE,  and  differs  from  the  mid-point  of  BF  by  the 
distance  AB ’ .  If  D  is  compared  with  the  images  B  and  F 
on  the  tilted  plate,  one  will,  therefore,  infer,  ignoring 
the  tilt,  a  direction  for  CD  that  is  in  error  by  the  angle 
(AB’)/(CD),  or  ,  nearly. 

Angular  errors  can  therefore  arise,  in  the  direc¬ 
tions  of  measured  images  inferred  from  the  star-images, 
equal  to  the  angular  error  of  the  base  of  the  plate  multi¬ 
plied  by  the  square  of  the  angular  semi-field  of  the 
plate.  The  semi-field  for  the  Goerz  camera  is  about  20°, 
or  one-third  of  a  radian;  consequently,  the  greatest  error 
in  an  inferred  direction  will  be  about  one-ninth  of  the 
angular  error  of  the  adopted  base.  If  the  zenith  is  adopted 
as  the  base,  then  an  error  of  less  than  1”  of  arc  will  result 
from  a  departure  of  the  plate  of  9”  from  horizontal.  This 
limit  of  9n  can  readily  be  met  by  careful  levelling,  causing 
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.the  greatest  errors  arising  from  uncertainty  of  the  base 
to  be  smaller  than  one  second  of  arc. 

5.  Corrections  for  Lens  Distortion  and  Astronomical 
Refraction.  The  camera  lens  distorts,  causing  images  to  be 
displaced  radially  away  from  or  towards  the  base  of  the  plate. 
For  any  but  a  very  poor  lens,  or  one  badly  tilted  in  its 
cell,  the  distortion  is  purely  radial;  and  the  effect  of 
atmospheric  refraction  is  likewise  (with  the  plates  '’based” 
at  the  zenith)  purely  radial  from  the  base.  Both  effects 
can  therefore  be  considered  at  once,  since  they  are  similar. 
Further,  it  is  easier  to  consider  effects  and  corrections  to 
the  standard  coordinates  X,  Y  than  to  the  plate  coordinates, 
measured  on  an  engine,  of  photographic  images.  The  standard 
coordinates  will  be  compared  eventually  with  the  plate 
coordinates,  and  therefore  it  is  logically  immaterial  whether 
one  corrects  the  plate  coordinates,  or  allows  in  the  standard 
coordinates,  for  distortion  and  refraction;  and  the  latter 
procedure  is  easier.  Atmospheric  refraction  moves  the  ob¬ 
served  position  of  a  star  towards  the  zenith  by  an  angle 

z  =  u  tan  z  (8) 

where  z  is  the  zenith  distance  and  * 


_  983b' 
460+ 1 


(9) 


where  Jl!_  is  the  barometer  reading  at  the  camera  in  inches, 
and  t  is  the  temperature  of  the  air  near  the  camera  in 
degrees  Fahrenheit.  These  formulae  represent  observed  re¬ 
fractions  correctly  within  a  second  of  arc  for  zenith  dis¬ 
tances  less  than  75°,  and  for  distances  under  20°  are  con¬ 
siderably  more  accurate.  At  Aberdeen,  u  is  about  58”  and 
this  value  can  be  taken  as  standard,  with  an  accuracy 
sufficient  for  present  purposes  (for  reasons  which  will  be 
obvious  later) .  It  is  easily  verified  that  within  20°  of 
the  zenith,  it  is  immaterial  whether  the  refracted  or  unre¬ 
fracted  tan  z  is  used  in  formulae  (8),  for  the  difference 
of  refraction  thereby  introduced  is  smaller  than  0”.01. 

Denote  by  r  the  distance  (X2  +  y2)1//2  of  the 
point  X,  Y  from  the  origin  on  the  standard  coordinate 
plane. 

*  Comstock,  Sidereal  Messenger.  April,  1890. 
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K 


0 


The  effect  of  refraction  is  to  diminish  z  by  ur/K,  and  thus 
to  diminish 


by 


r  =  K  tan  z 

dr  =  -  K(sec2  z)  ur/K 
=  -ur(l  +  r2/K2) . 


The  first  term,,  linear  in  r,  amounts  merely  to  a  change  of 
scale  and  can  be  ignored,  since  the  reduction  procedure  will 
automatically  allow  for  scale  factors.  We  are  left  with  the 
third-order  term 

dr  =  -ur3/K2  (10) 

where  u  must  be  exoressed  in  radians.  With  u  =  58”,  K  =  10, 

(10)  is 

dr  =  -0.000002812  r3  (11) 

This  amounts  only  to  2”. 7  at  the  edge  of  a  Goerz  plate,  and 
changes  only  0”.l  for  a  change  of  one  inch  in  the  barometer, 
or  of  17°  F.  It  is  for  the  preceding  reason  that  a  constant 
value  of  58”  may  be  adopted  for  u  at  the  Proving  Ground, 
with  sufficient  accuracy. 


The  effect  of  distortion  is  similar,  and  can  in¬ 
volve  only  odd  powers  of  r.  The  diagram  below  shows  the 
nodal  point  of  the  lens  at  0;  and  a  ray  with  an  angle  of 
incidence  of  z,  before  entering  the  lens,  and  with  an  angle 
of  incidence  zd_  after  traversing  the  lens  —  the  difference 
(z*  -  z)  being  the  angular  distortion.  It  is  clear  that  r_^ 
is  an  odd  function  of  _z,  and  hence  that  the  linear  distortion 
(r *  -  r)  must  be  an  odd  function  of  r  involving  only  odd 
powers  of  r.  . 
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For  the  same  reason  as  in  the  case  of  refraction,  the  terms 
linear  in  r  can  be  ignored,  being  absorbed  in  the  scale  fac¬ 
tor,  and  we  are  left  only  with  cubic  and  higher  odd  powers 
of  r,  thus: 

dr  =  A’  r^  +  B'r'3  +  ... 

In  combination  with  the  cubic  terms  in  refraction,  the  com¬ 
bined  effect  is  therefore 

dr  =  A”r^  +  B’r'5  +  ...  (12) 

The  constants  A”  and  B’  (B’  will  probably  be 
trivial)  must  be  found  once  and  for  all  empirically,  by 
the  method  to  be  described  shortly,  for  each  camera  lens. 
Although  A"  includes  refraction,  dependent  on  atmospheric 
conditions,  changes  in  the  part  of  A"  that  arises  from 
refraction  are,  as  we  have  seen,  minute  and  ignorable  and 
thus  A”  and  B»  may  be  regarded  as  constant  for  practical 
purposes.  The  correction  of  the  X  and  Y  of  a  star,  to 
allow  for  dr  given  by  (12),  is  carried  out  as  follows. 

Set  X  =  r  cos  0,  Y  =  r  sin  0,  where  0  is  an  azimuthal  .  ^ 

angle  (on  the  standard  coordinate  plane)  that  is  not  .  , 
altered  by  refraction  and  distortion.  Then 

dX  =  cos  0  dr;  dY  =  sin  9  dr 

and  since 

r2  _  xs  +  Ys 

one  has 


(dX)d  =  X  A"  (X2  +  Y2) 
(dY)d  =  Y  A”  (X2  +  Y2) 


1  +  W  (X2  +  Y2) 


1  +  f;  (X2  +  Y2). 


(XIII) 
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These  equations  enable  the  computer  readily  to  correct  the 
values  of  X  and  Y  already  computed  for  each  star,  thus: 


X«  =  X  +  (dX)d 
Y'  =  Y  +  (dY)d 


(XIV) 


where  the  Xf  and  the  YT  are  standard  coordinates,  with 
refraction  and  distortion  allowed  for. 


6.  Determination  of  the  Constants  A11  and  BT  (Distor¬ 
tion  and  Stellar  Refraction.  To  determine  A"  and  B*  for  a 
camera,  it  is  to  be  very  carefully  levelled  and  then  a  star- 
photograph  is  to  be  taken.  Compute  the  uncorrected  standard 
coordinates,  X  and  Y,  for  each  star.  A  ’'star”  can,  of  course, 
be  a  break  in  a  star-trail,  or  the  average  of  a  number  of 
breaks  in  the  same  trail.  By  the  method  later  to  be  described, 
find  by  least  squares  the  plate  constants  a,  b,  c,  d,  e,  f  in 
the  equations 


a  +  bx+cy  =  X 
d+ex+fy  =  Y 

where  x  and  y  are  plate  coordinates.  Then  find  the  constants 
A”  and  B*  from  the  residuals  by  least  squares,  using  the  ob¬ 
servational  equations 

C'X  +  A"X  (X2  +  Y2)  +  B’X  (X2  +  Y2)2=  a  +  bx  +  cy  -  X 

C»Y  +  A”Y  (X2  +  Y2)  +  B'Y  (X2  +  Y2) 2=  d  +  ex  +  fy  -  Y 

The  unknowns  are  C’,  A”,  and  where  C»,  although  it  must 

be  included  in  the  least  squares  solution,  is  of  no  interest 

when  found.  The  stars  should  cover  the  plate  as  fully  and  as 
uniformly  as  possible,  for  best  results,  and  the  present  deter¬ 
mination  need  be  carried  out  only  once  and  for  all,  for  each 
camera. 


7.  Determination  of  the  Plate  Constants.  The  equations 


a  +  bx  +  cy  =  X* 
d  +  ex  +  fy  =  Yr 


(XV) 


relate  the  measured  plate  coordinates,  x  and  y,  to  the  correc¬ 
ted  standard  coordinates  X'  and  Y* .  The  equations  (XV)  allow 
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for  scale  factor,  for  rotation  of  the  plate  in  the  measuring 
engine,  for  lack  of  perpendicularity  of  the  ways  of  the 
measuring  engine,  for  a  difference  between  the  scales  of  the 
x  and  y  screws,  and  for  arbitrary  zero-points  of  x  and  of  y. 
Three  stars  will  determine  the  plate  constants  a,  b,  c,  d,  e 
and  f-,  more  than  three  stars  will  overdeterraine  the  plate 
constants  and  a  least  squares  solution  should  be  made.  The 
solution  with  n  stars,  1,  2,  ...  n,  proceeds  as  follows. 

Let  x  be  the  mean  of  the  x*s,  y  be  the  mean  of  the  y’s,  and 

x»  =  x  -  x 

yt  =  y  _  y. 

Then  the  observational  equations  are 

(a  +  bx  +  cy)  +  bx*^  +  cy»^  =  X’^ 

(d  +  ex  +  fy)  +  ex’^  +  fy»^  =  Y’^ 

and  the  normal  equations  (the  unknowns  being  (a  +  bx  +  cy  = 
b,  c,  (d  +  ex  +  fy  =  d’ ) ,  e,  and  f)  are 

a  ’  b  c  = 

n  0  0  X» 

x»2  x«y»  x»X» 

y'2  y»x»  : 

with  similar  equations  in  Y*  for  dt,  e,  and  f.  To  solve 
these,  compute 

A  =  [x»2]  [y*2]  -  [x»y*]2 

bi  ={x’i  [y,s]  -  y*i  [x,y’] }  /A 

o±  =jy'±  [x’2]  "  x’i  [x»y»]J./A 

and  check  these  quantities  by  the  equations 


(XVI) 


at)  , 


(XVII) 
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1 - 1 

o' 

H* 

l _ 1 

it 

o 

o 

ii 

o 

1 _ 1 

b^Xf  -  1 

Vi}  >  0 

(XVIII) 

cixi]  =  0  j 

Vi]  =1-. 

Then 

a  =  (l/n)  X* i  -  bx  -  cy 

L . 

d  -  (X/n)  Y 1 ^  -  ex  -  fy 

b  =  biXf i 

e  =  [bp-! 

(XIX) 

c  =  M’i 

f  =  [c±Y»' 

and  when  only  three  stars  are  used,  equations  (XV)  may  be 
used  for  final  checks  and  will  be  satisfied  exactly.  With 
more  than  three  stars,  equations  (XIX)  must  be  repeated  to 
provide  a  check;  another  check  is  that  different  plates 
will  have  nearly  the  same  plate  constants  unless  the  cameras 
have  been  disturbed  between  exposures. 

When  there  are  only  three  stars,  the  preceding 
solution  holds,  but  it  is  quicker  merely  to  solve  the  six 
simultaneous  equations  (XV)  for  a,  b,  c’,  d,  e,  and  f  by 
successive  eliminations.  The  equations  group  into  two  sets 
of  three.  .. 


3.  Determining  the  Standard  Coordinates  of  the  Aerial 
Point.  One  has  merely  to  insert  in  equation  (XV)  the  plate 
coordinates  of  any  aerial  point  in  order  to  find  its  stand¬ 
ard  coordinates,  X*  and  Y*.  These  are  to  be  corrected  by 
applying  dX  and  DY  given  by  (XIII) ,  thus: 

X  =  xt  -  (dx)d  (XX) 

Y  =  Y>  -  (dY)d 

where  (XIII)  may  be  entered,  to  sufficient  accuracy,  with 
X*  and  YT  instead  of  X,  Y.  These  standard  coordinates  have 
next  to  be  corrected  for  the  difference  between  astronomical 
refraction,  and  the  smaller  refraction  that  affects  the  ray 
to  the  airplane  immersed  in  the  atmosphere. 


f 
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9.  The  Refraction  Affecting;  the  Aerial  Point.  If  the 
aerial  point  were  at  infinite  height,  its  refraction  would  be 
astronomical.  In  the  d^agram^  g-oo  _ _  I 


0  is  the  camera;  Z  is  the  zenith;  OPI  is  the  ray  from  the 
aerial  point;  P  is  the  aerial  point;  OQ  is  the  tangent  to 
the  ra}’-  at  0;  zQ  is  the  zenith  distance  of  the  ray  at  the 
camera;  z^  is  the  zenith  distance  of  the  straight  line  OP; 
Zee  is  the  zenith  distance  of  the  ray  OPI  extended  to 
infinity;  h  is  altitude  above  the  camera.  The  normal  dis¬ 
tance  from  00  to  the  ray  OPI  at  the  distance  s_  from  0  is  y. 
Now  if  z  is  the  variable  zenith  distance  along  the  ray  and 
if  u  is  the  index  of  refraction  of  air,  one  has 


or  since 


one  has 


dz  „ 
ds 


-  sin  z 


gii 

udh 


dh 

ds 


cos  z, 


cot  z  dz  = 


dy, 

U 


(21) 


\ 

) 
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whence 


sin  z  =  (  uQ/u  )  sin  zQ 


(22) 


where  u  is  the  index  of  refraction  of  the  air  near  the 
camera. 

Since  z  -  zQ  =  Az  is  small,  one  has 

Ll  r-'n  ' 

Az  =  tan  zn 

M'  O 


whence 


y  =  I  z  ds 
o 


=  tan  zQ  sec  z0  J  ff0--u  dh 


h 

J  K-U 


(23) 


and 


-  z  =  y/h  sec  z 


Jh  "  "o 


o 


tan  z  tt 
o  h 


h 

o 


U 


dh. 


(24) 


We  have 


w-  ~  1  +  (  u'o  “  -1)  P/Pc 


where  p  is  the  air  density,  so  that 


'  h 

zh  -  =  («*,  -  1)  tan  z0  i  '  J  1  -  (p/p0) 


dh.  (25) 


1+( H  0“D p/Po 


me  denominator  is  closely  enough  equal  to  unity,  since 
n  -  1  is  small;  and  thus,  very  nearly. 


-  z. 


=  ( 


1)  tan  z 


h 


1  - 


1 

h 


S  (p/p0) 


The  astronomical  refraction  is^  zx  -  z  ,  given  by  setting 
h  =  oo  m  (26):  ° 


z  -  zo  =  (  uQ  -  1)  tan  zQ  , 


which  agrees  with  (8)  if  u  =  u0  -  1  is  taken  as  53"  of  arc, 
or  as  0.0002812  radians.  Thus  the  difference  between  the 
astronomical  refraction  and  the  actual  is 


zoo“  zh  =  (1/1000)  tan  zQ  f(h) 


where 


h 

f(h)  =  0.2812  i  S  (p/pQ)  dn. 


The  value  of  f(h)  up  to  130,000  feet  has  been  com 
puted  from  the  annual  means  of  the  observed  values  of  p/p  , 
obtained  by  sounding  balloons,  and  published  by  Humphreys?* 


h 

(ft) 

f(h) 

h 

(ft) 

f(h) 

h 

(ft) 

f(h) 

0 

.281 

45000 

.144 

90000 

.083 

5000 

.258' 

50000 

.135 

95000 

.079 

10000 

.239 

55000 

.126 

100000 

.075 

15000 

.222 

60000 

.118 

105000 

.072 

20000 

.206 

65000 

.110 

110000 

.069 

25000 

.191 

70000 

.104 

115000 

.066 

30000 

.178 

75000 

.098 

120000 

.063 

35000 

.166 

80000 

.093 

125000 

.061 

40000 

.155 

85000 

.087 

130000 

.059 

*  W.  J.  Humphreys,  Physics  of  the  Air.  Franklin  Institute,  1920. 
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The  above  values  apply  to  a  barometer  reading  of  30  inches, 
and  a  temperature  of  43°  F,  at  the  camera.  In  accordance 
with  equation  (9) ,  they  should  be  lowered  by  1%  for  each  5°  F 
by  which  the  temperature  exceeds  43°  F,  and  increased  by  1% 
for  each  0.3  inches  by  which  the  pressure  exceeds  30  inches. 
The  altimeter  will  furnish  values  of  h  accurate  enough  for 
finding  f(h)  by  means  of  the  preceding  table  and  corrections. 

By  the  argument  of  section  5,  we  allow  for  the 
difference  of  refractions  by  applying  the  correction  dr, 
where 

1000  dr  =  -K(l  +  r2/K2)  r  f(h)/K 
or  by  applying  the  corrections  (dx).„,  (dY)^,  where 

1000  (dX)r  =  f (h) 


1000  (dY)r  =  f (h) 


1  +  (X2  +  Y2)/K2 


1  +  (X2  +  Y2)/K£ 


] 


(XXVIII) 


to  the  X  and  Y  of  the  aerial  point,  in  order  to  find  X1',  Yn : 


X”  =  X  -  (dX)p 
Y"  =  Y  -  (dY) r* 


(XXIX) 


10.  The  Direction  Ratios  to  Aerial  Points.  With  respect 
to  the  set  of  axes  X,  Y,  Z,  defined  in  section  3  (X  to  the 
east,  Y  to  the  north,  and  Z  to  the  zenith) ,  having  an  origin 
at  either  camera:  The  geometrical  coordinates  of  the  inter¬ 
section,  of  the  standard  coordinate  plane  with  the  straight 
line  from  the  camera  to  the  aerial  point,  are  simply  X,!,  Y” ,  K. 
Hence  the  direction  ratios  of  the  line  from  the  camera  to  the 
aerial  point  are  merely 


.  X”  :  Y"  :  K.  (XXX) 

It  will  be  henceforth  necessary  to  distinguish  between 
the  two  cameras,  which  will  be  denoted  as  camera  1  and  camera  2. 
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The  ratios  for  camera  1  will  be  denoted  by  X'']_  :  Y11!  :  K 
and  those  for  camera  2  by  X,r2  :  Yf,2  :  K.  It  will  be  neces¬ 
sary  to  transform  these  direction  ratios  to  a  standard 
terrestrial  system  of  reference.  Because  of  the  curvature 
of  the  earth,  the  X,  Y,  Z  systems  of  the  two  cameras  will 
not  be  parallel  to  each  other.  If  the  standard  system  of 
reference  were  chosen  to  coincide  with  the  X,  Y,  Z  system 
of  one  of  the  cameras,  some  computation  would  be  saved; 
but  then  the  direction  of  gravity  over  the  bomb-trajectory 
would  not  coincide  with  the  standard  z-direetion.  We  dis¬ 
cuss  the  transformations  as  though  the  standard  reference 
system  differed  from  the  X,  Y,  Z  systems  of  both  cameras; 
if  it  is  identical  with  either  of  the  latter  then  the 
following  results  and  procedure  will  still  hold. 

Let  Pp  denote  the  origin,  at  latitude  cp0,  of  a 
standard  terrestrial  reference  system,  §  ,  q  ,  C .  The 
latitudes  of  cameras  1  and  2  are  denoted  by  cp]_  and  ©2,  and 
their  longitudes  measured  to  the  west  from  PD  are  denoted 
by  ai  and  a  .  XG,  Yg,  ZQ  are  a  set  of  axes  through  P0  as 
origin,  sucli  that  ZQ  is  vertical,  Y  is  to  the  north,  and 
X0  is  to  the  east.  Xi,  Y]_,  Zi  and  X2,  Y2,  Z2  are  similar 
sets  of  axes  through  cameras  1  and  2  as  origins;  these 
three  sets  of  axes  are  of  course  not  parallel  to  each  other. 
The  standard  reference  system  differs  from  the  Xo>  Y0,  ZQ 
system  through  a  rotation  about  the  ZQ-axis,  which  coincides 
with  the  C  -axis.  Then  -axis  has  an  azimuth  at  PQ  of  A  , 
measured  from  north  to  east;  the  E  -axis  has  an  azimuth  at  PG 
of  (  A  +  90°) .  We  also  make  use  in  this  discussion  of  the 
hour-angle  declination  axes  x’Q,  y*0,  z»D;  with  P0  as  origin, 
with  x»0  lying  in  the  equator  ana  meridan,  y»0  pointing  to 
the  east,  and  z»0  pointing  to  the  north  celestial  pole. 
x|l>  Y’l,  and  xt2>  Y’2>  z' 2  are  similar  axes  with  ori¬ 
gins  at  cameras  1  and  2,  respectively. 

It  is  necessary  to  find  the  cosines  of  the  angles 
between  the  E  ,  q  ,  £  axes  and  the  Xi  ,  Yi,  Z±  axes;  and 
between  the  E  ,  q  ,  £  axes  and  the  X2,  Y2,  Z2  axes.  The 
cosines  of  the  angles  between  the  xQ  and  the  £  systems  are 


E 


n 


c 


all  =  cos  ^ 

al2  =  “sin  ^ 

a  =0 

13 

a21  =  sin  ^ 

a^  =  cos  X 

a23  =  0 

a3!  =  0 

a32  =  0 

a^^  =  1 
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The  cosines  relating  the  system  to  the  x'  system  are 


bn  =  o 

b12  =  1 

b13  =  0 

b21  =  -sln  ‘Po 

b22  =  0 

b23  =  cos  “’o 

b31  =  cos  cpQ 

b32  =  0 

b33  =  sln  1>o 

the  cosines  relating  the  x^  system  to  the  x»n  system  are 


cl.l  =  cos  “t 

c12  =  sdn  ai 

. 

o 

ii 

H 

O 

°21  =  -sln“, 

cno  =  cos  a 

22  i 

°23  =  0 

C31  =  ° 

o 

?o 

II 

o 

c33  =  1 

and  the  scheme  relating  the  systems  x’^  and  X-j_  is 

Yn  Z, 


dll  "  0 


d21  “  1 


=  -sin  (p^_  =  cos 


d22  0 


d23  “  0 


d32  =  cos  ^1  d33  =  sdn  ^1 


0 


£ 

n 

K 

where 


gll 

g12 

g13 

g21 

g22 

s23 

g31 

g32 

g33 

fik  dkj  . 


The  preceding  computations  are  easily  performed  if  actual 
numbers  are  inserted  in  the  cells,  and  lead  to  the  numeri¬ 
cal  array,  g-;  relating  the  Xi,  Yj,  Zi  system  to  the 
standard  system.  Then  the  numbers  and  headings  should  be 


X"i 

Y"l 

+ 

£ 

gll 

g12 

Kg13 

h 

g21 

g22 

Kg23 

C 

g31 

S32 

Kg33 

written  down  on  a  heavy  card,  headed  ’'Camera  1”. 

In  precisely  the  same  way,  replacing  camera  1 
by  camera  2,  and  thus  cpi  by  cp2  and  o.  by  oc  ,  one  obtains 
the  scheme,  which  should  be  written  &own  numerically  on  a 
heavy  card,  headed  "Camera  2". 


X"2 

Y"2 

+ 

£ 

hll 

h12 

Kh13 

h 

h2l 

h22 

^23 

4 

h31 

h32 

Kh33 

-24- 


The  transformation  of  the  direction  ratios 
:  Y"]_  K  to  the  system  of  standard  terrestrial  axes 
5,7],  £  is  merely 

A1  =  X"l  gll  +  y,'l  g12  +  (Kg13} 

BX  =  X"i  g2i  +  y"i  g22  +  (Kg23) 

°X  =  X”l  S31  +  Y"l  S3 2  +  (K«33) 


where  A]_  :  B]_  :  Ci  are  the  direction  ratios,  of  the 
straight  line  from  camera  1  to  the  aerial  point,  expressed 
in  the  £  ,  r  >  £  system.  Likewise,  the  direction  ratios 

X,! 2  '  Y’!2  :  K  from  camera  2  are  transformed  to  the  ratios 
A2  :  B2  :  C2  in  the  ^  ,  rj  ,  (  system  by  the  equations 


A2  =  X”2  hll  +  Y,!2  hl2  +  ^Kh13} 


B0  =  X”2  h21  +  Y”2  h22  +  (Kh97) 


23' 


C2  Xn2  h31  +  Y"z  h32  +  ^Kh33^ 


Using  the  coefficients  written  on  the  cards,  a  good  com¬ 
puter  should  be  able  to  transform  a  single'  set  of  direc¬ 
tion  ratios  X”  :  Y”  :  K  in  little  more  than  one  minute; 
the  coefficients,  which  depend  only  on  the  camera  posi¬ 
tions  and  on  the  standard  reference  system  that  is  used, 
are  mere  constants. 

If  the  origin  P0  of  the  reference  system  5  ,  , 

£  coincides  with  either  camera,  and  if  the  q  -axis  runs 
exactly  north  from  that  camera,  then  no  transformation  is 
necessary  of  the  ratios  X>'  :  Y”  :  K  obtained  from  that 
camera,  and  one  has  for  that  camera  A  =  X”  :  B  =  Y”  :  C  =  K. 
Thus  some  saving  of  computation- time  can  result  from  such  a 
choice  of  a  standard  reference  system. 


11*  The  Coordinates  of  the  Aerial  Point.  Denote  the 
£  >  7]  >  £  coordinates  of  camera  1  by  &1,  t^,  c^  and  those 


(XXXI) 


(XXXII) 
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of  camera  2  by  £2,  bg,  C2*  Denote  the  direction  ratios  of 
the  rays  from  the  two  cameras  to  the  aerial  point,  in  the 
x,  y,  z  system,  by  Aq  :  Bq  :  Cq  and  A2  :  %2  ;  ^2  as  in 
section  10.  Then  the  coordinates  x,  v,  z  of  the  aerial 
point  are  given  by  the  equations 


£  =  a-,  +  An  p  =  a0  +  A0  m 


T)  =  +  Bi  ?  =  +  B2  m 


C  =  C1  +  Cj  p  =  c2  +  C2  m 


(XXXIII) 


where  p  and  m  are  found  by  solving  the  observational  equations, 
by  least  squares. 


p  m  = 


Compute 

F  =  A\  +  B21  +  C2^ 

G  =  a2z.+  b2z  +  c22; 

H  —  A -^2  +  ®1®2  ^1^2’ 

I  -  (ag- a-^)  A-^  +  (h2— b^)  B1  1 

J  =  (a2_a1)  A2  +  (b2-bi)  B2  +  (c2-c1)  C2j 


(34) 


(XXXV) 
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Then 


p  =  (GI  -  HJ)/K 
m  =  (HI  -  FJ)/K 


(XXXVI) 


and  these  computations  should  be  checked  by  computing 
residuals  of  the  three  observational  equations,  and 
seeing  that  the  sum  of  their  products  with  the  coeffi¬ 
cients  of  should  vanish,  and  that  the  sum  of  their 
products  with  the  coefficients  of  m  should  also  vanish. 
This  checks  the  least  squares  solution.  Then  from 
(XXXIII),  one  finds  two  values  of.^'5  ,  of  r]  ,  and.  of  £  ; 
the  agreement  indicates  the  accuracy  of  the  whole  work, 
and  for  final  values  the  averages  of  the  two  estimates 
of  each  quantity,  £  ,  rj  ,  and  £  ,  should  be  taken. 

This  completes  the  description  of  the  method 
of  finding  the  position  of  an  aerial  point,  from  photo¬ 
graphs  taken  from  two  vertical  cameras,  by  the  use  of 
star-images.  We  can  call  the  preceding  method  the  "long 
method". 


12.  Summary  of  the  "Long  Method";  Directions  for 
Computers;'  "Estimated  Computation-Times .  The  preceding 
exposition  has  been  long,  but  the  method  is  fairly 
rapid.  It  may  be  summarized  as  follows.  At  the  outset, 
a  study  should  be  made  as  described  in  section  6  to 
find  the  distortion  and  third-order  refraction  terms  for 
the  two  cameras.  At  the  outset,  also,  when  the  piers 
are  set  out,  one  should  prepare  the  two  cards,  each  with 
nine  constants,  described  in  section  10.  Then  the  pro¬ 
cedure  v/ith  two  plates,  one  from  each  camera,  in  a  parti¬ 
cular  range-bombing  program  is  as  follows:  Measure*  the 
plate  coordinates,  (x,  y) ,  of  each  comparison  star  and 
aerial  point  on  both  plates.  While  this  is  being  done, 
compute  the  direction  cosines  Lp,  M0,  N0  of  the  comparison 
stars  by  equations  (i)  and  (ill)  --  unless  they  are  al¬ 
ready  available  from  the  Computations  of  a  previous  range 
bombing.  Correct  these  for  nutation,  etc.,  by  equations 
(IV)  and  compute  their  standard  coordinates  for  each 


-27- 


camera  by  eouations  (VI)  and  (VII) .  Apply  the  checks 
described  in  the  text,  at  all  stages.  Apply  the  correc¬ 
tions  given  by  equations  (XIII)  to  all  standard  coordi¬ 
nates,  by  equations  (XIV).  Determine  the  plate  constants 
from  the  comparison  stars,  using  equations  (XVI )  through 
(XIX) .  Determine  the  standard  coordinates  of  the  aerial 
points  by  equations  (XV),  correct  them  by  equations  (XX) 
and  (XXIX) ,  and  compute  their  direction  ratios  by  equations 
(XXXI)  and  (XXXII),  using  the  two  cards.  Finally  compute 
the  space  coordinates  of  the  aerial  points  by  using 
equations  (XXXV),  (XXXVI),  and  (XXXIII) .  The  writer  esti¬ 
mates,  very  roughly,  a  computation-time  of  from  four  to 
eight  hours  for  a  single  computer  when  there  are  three 
comparison  stars,  and  three  aerial  points  to  locate. 

With  a  suitable  arrangement  of  the  work,  three  computers 
should  accomplish  it  in  about  two  to  four  hours.  These 
times  are  for  experienced  computers;  mediocre  ones  would 
probably  take  longer. 


13.  Rapid  Method  Without  Using  Stars.  If  the  cameras 
are  very  rigidly  mounted  on  well-settled  piers  (protected 
further  from  diurnal  or  other  rapid  temperature  changes, 
and  not  moving  with  the  tides) ,  it  may  prove  possible  to 
avoid  the  employment  of  star-images  on  some  of  the  plates. 
For  this  "short  method"  to  be  possible,  fiducial  marks  must 
be  impressed  on  the  plates  by  marking  devices  that  are 
absolutely  fixed  with  respect  to  the  lenses  and  piers. 
Further,  the  lens  of  each  camera  must  be  absolutely  free 
from  rattle  in  its  cell,  and  the  camera  itself  must  be  of 
the  firmest  and  most  rigid  construction.  These  require¬ 
ments  are  unnecessary  when  star  images  are  used,  as  in 
the  "long  method". 

Using  stars  as  already  described,  it  is  necessary 
to  find  the  plate  constants  a,  b,  c,  d,  e,  and  f  (as  ex¬ 
plained  in  paragraph  7)  for  some  plate  which  we  shall  call 
the  "standard"  plate,  which  need  not  photograph  the  aerial 
point.  Coordinates  of  points  on  the  standard  plate  will 
be  written  thus:  x,  y.  Then  the  plate  constants  relate 
x  and  y  to  directions  in  space.  Suppose  no?/  that  we 
measure  the  coordinates  x,  y  for  a  number  of  fiducial 
marks  impressed  on  the  standard  plate  by  the  camera,  and 
then  take  another  plate,  (the  "observing"  plate),  on  which 
the  same  marks  are  impressed,  v/ithout  disturbing  the  camera 
in  any  way  betv/een  exposures.  We  can  denote  coordinates 
on  the  observing  plate  by  (x!  ,  y*)  and-*  we  can  measure  the 

*  The  plate  coordinates  (x1,  y’ )  are  not  the  space  coordinates 
(x1,  y1,  z')  of  Section  3. 


coordinates  (x’  ,  y’)  of  the  fiducial  marks  and  of  "the 
aerial  points  on  the  observing  plate.  Ho  stars  need 
appear  on  the  observing  plate.  Then  it  is  obvious 
that  if  there  are  three  or  more  fiducial  marks  on 
both  plates,  we  can  find  from  the  coordinates  x,  y  and 
x» ,  y’  of  the  fiducial  marks,  the  six  coefficients  in 
two  linear  equations,  like  the  equations  (XV) ,  which 
relate  the  coordinates  x»  y*  to  the  coordinates  x,  y. 

Such  relations  will  automatically  allow  for  changes 
of  scale  in  the  measuring  engine,  for  lack  of  perpen¬ 
dicularity  of  the  ways,  for  rotations  of  the  plates  in 
the  measuring  engine,  and  for  differences  between  the 
scales  of  the  two  screws  of  the  engine.  They  will  also 
automatically  allow  for  thermal  expansions  of  the  plates, 
and  of  the  camera.  Thus  from  the.x’ ,  y’  of  an  aerial 
point  on  the  observing  plate,  we  can  compute  its  x  and  y 
by  the  linear  relations'  so  found,  and  thence  deduce  the 
directions  X”  :  Y”  :  K  of  the  line  to  the  aerial  point, 
and  thence  finally  obtain  the  space  coordinates  B,  ,  q  ,  £ 
of  the  aerial  point,  just  as  though  the  aerial  point  had 
been  photographed  on  the  standard  plate  itself.  Care 
must  be  taken  to  allow  for  differences  between  the  refrac¬ 
tions  influencing  the  two  plates,  if  they  are  taken  on 
different  days;  but  the  procedure  is  straightforward  and 
it  is  left  to  the  reader  to  develop  the  formulae  appli¬ 
cable  to  the  case  of  three  or  more  fiducial  marks. 
Unfortunately,  with  the  Goerz  cameras  of  the  Aberdeen 
Proving  Ground,  there  are  only  two  fiducial  marks  on 
each  camera.  We  shall  discuss  this  case  with  care. 

Let  us  inquire  as  to  what  differences  there  can 
be  between  the  standard  plate  and  the  observing  plate, 
apart  from  changes  of  refraction,  when  the  camera  has  not 
been  disturbed.  The  xL? s  and  xL' s  can  differ  from  the 
x’s  and  y’ s  through  changes  in  temperature,  either  at  the 
times  of  measurement  or  at  the  times  of  exposure,  or  both. 
Changes  of  temperature,  whether  affecting  the  plates  through 
their  glass,  or  through  swelling  or  contracting  of  the 
camera,  or  through  changes  in  the  pitches  of  the  measuring 
screws,  can  at  most  amount  to  a  simple  change  in  scale. 

Thus  we  should  expect  the  two  plates  to  differ  in  scale. 

If  the  plates  are' measured  on  different  engines,  then  the 
pitches  of  the  two  x-screws  can  differ  in  ratio  from  the 
pi tones  of  the  two  y-screws.  This  must  be  avoided,  when 
only  two  fiducial  marks  are  available,  by  using  the  same 
engine  for  measuring  both  plates.  Further,  we  must  take 
care  in  the  present  ’’short”  method  not  to  use  an  engine 
which  has  only  one  screw,  and  y/hich  depends  for  measure¬ 
ments  of  x  and  y  upon  a  rotation  of  the  plate  through  90° 
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between  measuring  the  x* s  and  measuring  the  y’s.  It  is 
essential  to  use  an  engine  with  two  screws,  very  nearly 
at  right  angles.  With  a  single-screw  engine  the  angle 
between  the  axes  of  measurement  will  not  only  never  be 
90°,  but  it  will  vary  from  plate  to  plate  according  to 
the  accuracy  with  which  the  x  and  y  dispositions  can  be 
reproduced.  While  that  has  no  evil -effect  when  there 
are  three  fiducial  marks,  it  introduces  insuperable 
difficulties  when  there  are  but  two  marks.  With  a  two- 


screw  engine,  the  angle  between  the  ways  may  not  be  90° 
exactly,  and  in  general  it  will  not  be  90°;  but  whatever 
the  angle  is,  it  will  remain  constant  from  plate  to 
plate  unless  the  engine  is  subjected  to  abuse.  Moreover, 
the  ratios  between  the  scales  of  the  two  screws  will  re¬ 


main  constant  from  plate  to  plate,  since  the  screws  are 
always  made  of  the  same  metal.  Another  source  of  differ¬ 
ence  between  (x1 ,  y1)  and  (x,  y)  arises  from  the  inevit¬ 
able  small  differences  between  the  orientations  of  the 


plates  in  the  measuring  engine  when  the  plates  are 
measured.  These  cannot  be  entirely  avoided.  Finally, 
the  xM  s  and  y1 1 s  can  differ  from  the  x's  and  v’s  through 
the  "two  systems1  having  different  zero-points.^  We  are 
left,  therefore,  when  all  precautions  have  been  taken, 
with  the  effects  of  differences  of  scale,  of  different 
orientations  in  the  measuring  engine,  and  of  different 
zero-points.  We  discuss  the  relation  between  x1 ,  y» 
and  x,  y  in  the  following  section,  and  show  how  the 
effects  of  the  preceding  differences  can  be  automatically 
eliminated. 


M*  Short  Method  -  The  Relation  of  the  Standard  Plate 
to.  the  Observing  Plate.  In  the  diagram  below,  we  have, 
drawn  the  x1  and  y1  axes,  and  also  a  set  of  orthogonal 
axes  (X«»,  y111)  such  that  x”1  is  parallel  to  x1  .  We  make 
the  origins  coincide  here  and  in  the  following  transforma¬ 
tions,  since  zero-point  constants  are  easily  inserted 
later  on.  One  unit  of  x1  will  be  bs  units  of  x1’1  or  y" » 
and  the  angle  a  measures  the  lack  of  perpendicularity  of* 
the  ways  of  the  engine,  x1  and  yt .  A  unit  of  the  y1  screw 
is  cs  units  of  x» 1  or  y« ». 

/T’ 

yu  /  ' 

/  X1 

#/ - ^(xsy') 
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Then 


x»  ’  =  b  x»  +  c_y’  sin* 

S  5 


y” 1  =  c0y‘  cos* 

b 

The  total  change  of  scale  and  the  effect  of  the  different 
orientations  of  the  plates  in  the  engine,  can  be  repre¬ 
sented  by  a  transformation  from  (x,Tl,  y,M)  to  (x11,  y")  where 
the  latter  coordinates  differ  from  the  former  through  a 
rotation  oo  and  through  a  scale  factor  s..  Thus 


x"  =  sxJTI  cos  co  -  sy,T?  sin  o> 


yti  =  sx"*  sin  co  +  sy,!  f  cos  co. 


Finally,  the  measurement  of  the  standard  plate  in  the 
engine  corresponds  to  the  diagram 


where  the  relation  of  (x".  y")  to  (x,  y)  is  the  same  as 
the  relation  of  (x« ’ ,  yn»;  to  (x> ,  y»).  Thus 


x  =  _I_  (Xn  _  yu  tan  *) 
s 


_ Z!! _ 

cg  cos  a 


I 
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It  will  be  noticed  that  we  have  assumed  the  same  relative 
pitches  of  the  x  and  y  screws  for  the  two  plates,  and  the 
same  error  of  the  ways  «;  and  that  changes  of  scale  from 
one  plate  to  the  other  have  been  absorbed  in  the  constant 
s.  Combining  the  preceding  linear  equations,  and  adding 
zero-point  constants,  one  finds  that  in  general 


c  s 

X  =  r  +  s(c os  a)  -  sin  co  tan  a)  x1 — sin  a>  sec  u  .  y’ 


sin  co  sec  a 


.  x’  +  s(cos  co  +  sin  co  tan  k)  y* 


which  can  be  written 


x  =  r  +  Ax1  +  Q  y1 


+A  tan  co  tan  a 


y  =  n  -  QX«  +  Ay’  + 


+  Q  (l— ~ )  -xl  +A  tan  co  tan  a  .  y» 


where  A  has  been  written  for  s  cos  co  and  Q  has  been 
written  for  -s  sin  co  sec  a. 


Since  x  is  a  very  small  angle,  and  co  is  small, 
while  the  scale  factors  s,  bs,  and  cs  are  all  very  close 
to  unity,  it  follows  that  A  is  nearly  unity  while  Q  is 
approximately  equal  to  -sin  co.  The  ratios  of  scales  bs/cs 
is  very  close  to  unity  for  any  well-made  engine;  certainly 
not  differing  from  unity  by  more  than  one  part  in  10,000. 

The  angle  a  for  a  well-made  engine  should  be  less  than  one 
minute,  so  tan  a  is  less  than  0.0003.  If  (by  careful  aline- 
ment  of  the  observing  plate  in  the  measuring  engine  so  that 
x,  y  and  x’ ,  y*  are  nearly  the  same  for  both  fiducial  ooints) 
we  can  keen  co  less  than  1°,  then  the  coefficients  of  x»  and  y’ 
in  the  square  bracketed  second-order  terms  are  all  smaller 
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than  0.000005  and  are  ignorable.  We  are  then  left  with 

f  +  A  x*  +  q  y».  =  x 

n  -  c  x>  .+  a  y1  =  y 

with  an  accuracy  sufficient  to  yield  directions  that  will 
he  correct  to  better  than  a  second  of  arc.  Here  ,  A  , 

S  ,  n  are  unknown  constants  to  be  determined  from  the 
measured  coordinates  (x,  y)  and  (x> ,  y»)  of  the  fiducial 
marks.  It  is  easy  to  aline  both  plates,  in  the  measuring 
engine,  similarly  to  within  one  degree  of  angle.  We  may 
call  f  ,  A  ,  Q  ,  FI  the  "observing  constants". 


15 .  Short  Method  -  Finding  the  Observing  Constants. 
It  is  readily  found  that 


A 

= 

(* 

'i  -  x,2)(>i  -  x2)  +  fr'i  -  y'2)(n  -  y2}_ 

/E 

i 

r— 

— i 

Q 

l 

(y 

’i  -  y'2)(xx  -  x2)  -  (n  -  -  x'z> 

1  /E 

! 

r 

=  n  - 

-  Cy»x 

n 

=  n  + 

2x'-]_  -  Ay 

■ 

where 

E  =  (x»i  -  x<2)2  +  (y»i  -  y»2)2. 

The  observing  constants,  E  ,  A  ,  Q  , H  may  be  checked 
by  the  equations  (XXXVII) .  Here  (x’l,  y’l)  and  x’2*  y’2) 
are  the  measured  coordinates  of  the  fiducial  marks  1  and  2, 
respectively,  on  the  observing  plate;  and  (xi,  yi)  and 
(x2>  Y2)  are  the  measured  coordinates  of  the  fiducial  marks 
on  the  standard  plate.  As  has  been  pointed  out,  both  plates 
must  be  measured  in  the  same  two-screw  measuring  engine. 


(XXXVII) 


t  ■ 

(XXXVIII) 
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after  having  been  inserted  in  nearly  parallel  orientations. 
The  observing  constants  are,  of  course,  relative  to  the 
particular  observing  plate,  and  particular  standard  plate. 
It  is  entirely  proper  to  use,  as  a  standard  plate,  the 
average  of  two  or  more  star-plates  taken  at  nearly  the  same 
time;  then  one  has  to  use  (later) the  average  of  the  air 
temperatures  and  pressures  when  correcting  for  refraction 
on  the  standard  plate. 


16.  Short  Method  -  The  Direction  Ratios  to  Aerial 
Points.  Let  the  plate  coordinates  of  an  aerial  point  be 
(xT , y* )  as  measured  on  the  observing  plate.  Insert  these 
values  in  eauations  ( XXXVII)  and  thereby  obtain  (x,  y) . 

Now  insert  (x,  y)  in  the  left-hand  members  of  equations 
(XV) ,  involving  the  plate  constants  of  the  standard  plate, 
and  thereby  find  the  standard  coordinates  (X1 ,  Y’)>  of 
the  aerial  point.  The  original  values  (x’ ,  y»)  contained 
the  effects  of  refraction  of  the  aerial  point,  and  lens 
distortion.  The  refraction  in  question  is  that  appropriate 
to  the  date  of  the  observing  plate.  It  is  characteristic 
of  the  transformations  (XV)  that  they  automatically  take 
out  the  first  order  astronomical  refraction  for  the  date 
of  the  standard  plate.  Therefore,  the  standard  coordinates 
(X* ,  Y')  contain  the  refraction  of  the  aerial  point,  dis¬ 
tortion,  and  minus  the  first  order  astronomical  refraction 
of  the  standard  plate.  If  we  subtract  from  X*  and  Y*  the 
small  quantities  (dX)d  and  (dY)d  given  by  (XIII)  entering 
these  equations  with  (XT ,  Y’) ,  we  are  left  with  standard 
coordinates  free  from  distortion,  and  affected  by  minus 
the  full  astronomical  refraction  of  the  standard  plate, 
and  plus  the  refraction  of  the  aerial  point  on  the  date 
of  the  observing  plate.  If  we  now  apply  the  difference 
(astronomical  refraction  on  date  of  standard  plate)  minus 
(astronomical  refraction  on  date  of  observing  plate)  and 
apply  further  the  difference  (astronomical  refraction) 
minus  (refraction  to  height  h)  appropriate  to  the  date  of 
the  observing  plate,  we  shall  be  left  with  standard  coor¬ 
dinates  (X,!,  Y”)  fully  corrected  for  distortion  and  refrac¬ 
tion  and  comparable  in  all  respects  with  the  standard 
coordinates  (X”,  Yn)  furnished  by  equations  (XXIX)  of  the 
long  method.  The  difference  between  the  astronomical  refrac 
tions  on  the  two  dates  is  obtained  readily  from  equations 
(8)  and  (9) ,  and  thus  we  have  for  any  aerial  point 
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(XXXIX) 


X"  =  X'  -  (dx)d  +  (uob  -  ust)  X>  -  (dX)p 


Y'>  =  Y>  -  (dy)d  +  (uob  -  ust)  Y'  -  (dY)p 


where 


0.004766  b* 

460  +  t 


and  where  u0b  refers  to  the  observing  plate,  ust  to  the 
standard  plate.  Here  b_[_  is  the  air  pressure  at  the 
camera  in  inches  of  mercury,  and  t  is  the  air  temperature 
in  degrees  Fahrenheit. 

Thus  the  short  procedure  goes  as  follows:  In¬ 
sert  the  plate  coordinates  (xf  ,  y»)of  the  aerial  point, 
as  measured  on  the  observing  plate, 'into  equations  (XXXVII) 
and  thus  obtain  (x,  y) .  Insert  (x,  y)  into  the  left-hand 
members  of'  equations  (XV)  appropriate  to  the  observing 
plate,  and  thus  find  (X* ,  Y ’ ) .  Correct  (Xr,  Yf)  by 
equations  (XXXIX)  to  find  (X”,  Y") .  In  equations  (XXXIX) , 
find  (dX)d  and  (dY) a  from  equations  (XIII) ,  ignoring  the 
difference  between  (X,  Y)  and  (Xf ,  Y')>  ana  similarly 
find  (dx)r  and  (dY)r  from  equations  ( XXXVIII)  and  the 
associated  table. 


17.  Short  Method  -  The  Coordinates  of  the  Aerial 
Point.  The  quantities  (X",  Y")  just  obtained  are  precisely 
comparable  to  the  same  quantities  in  the  long  method,  and 
from  there  on  the  reduction  goes  just  as  in  the  long 
method,  described  in  sections  10  and  11. 


18.  Summary  of  the  ’’Short”  Method;  Directions  for 
Computers;  Estimated  Comoutation-Times.  Standard  plates 
containing  star-images  must  have  been  previously  obtained 
from  both  cameras;  then  the  cameras,  without  being  in  any 
way  disturbed,  photograph  the  aerial  points  on  observing 
plates.  The  comparison  star-images  and  fiducial  marks 
must  be  measured  on  the  standard  plates;  and  the  plate 
constants  must  be  obtained  as  in  the  "long"  method  for 
the  standard  plates.  The  fiducial  marks  and  the  images 
of  the  aerial  points  are  measured  on  the  observing  plates, 
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by  the  use  of  the  sane  engine  that  was  used  to  measure 
the  standard  plates.  The  standard  plates  and  observing 
plates  must  be  inserted  in  the  measuring  engine  (which 
must  be  of  the  two-screw  type)  approximately  parallel 
to  each  other.  Find  the  observing  constants  by  equations 
( XXXVI II )  ,  and  then  by  equations  ( XXXVI I)  ,  compute  stand- 
ard  plate  coordinates  (x,  y)  for  each  aerial  point. 

Compute  standard  coordinates  (X?  ,  Y')  by  equations  (XV) .  \ 

Correct  these  by  equations  (XXXIX) ,  obtaining  the  correc- 
ted  standard  coordinates  (X” ,  Y”)  of  the  aerial  points. 

By  formulae  (XXXI)  and  ( XXXII ),  compute  the  direction 
ratios,  using  the  two  cards.  Finally  compute  the  space 
coordinates  £  ,  p  ,  £  of  the  aerial  ooints  by  using 

equations  (XXXV) ,  (XXXVI),  and  (XXXIII) .  The  writer 
estimates  very  roughly  a  computation-time  of  from  three 
to  four  hours  for  a  single  computer  when  there  are  three 
aerial  points  to  locate.  With  suitable  division  of  the 
work,  three  computers  should  accomplish  it  in  about 
one  to  two  hours.  These  times  are  for  experienced  com¬ 
puters;  mediocre  ones  or  beginners  would  take  longer, 
in  connection  with  these  estimates,  the  standard  plates 
are  supposed  already  to  have  been  reduced. 


T.  E.  Sterne, 

Captain,  Ordnance  Dept. 
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